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Abstract 

We propose a phenomenological model of the Dirac neutrino mass matrix based on the Fridberg- 
Lee neutrino mass model at a special point. In this case, the Fridberg-Lee model reduces to the 
Democratic mass matrix with the S3 permutation family symmetry. The Democratic mass matrix 
has an experimentally unfavored degenerate mass spectrum on the base of tribimaximal mixing 
matrix. We rescue the model to find a nondegenerate mass spectrum by adding the breaking mass 
term as preserving the twisted Fridberg-Lee symmetry. The tribimaximal mixing matrix can be 
also realized. Exact tribimaximal mixing leads to ^13 = 0. However, the results from Daya Bay and 
RENO experiments have established a nonzero value for 613. Keeping the leading behavior of U 
as tribimaximal, we use Broken Democratic neutrino mass model. We characterize a perturbation 
mass matrix which is responsible for a nonzero 0\3 along with CP violation, besides the solar 
neutrino mass splitting has been resulted from it. We consider this work in two stages: In the 
hrst stage, we obtain the perturbation mass matrix with real components which breaks softly the 
/i — r symmetry and this leads to a nonzero value for 613. In the second stage, we extend the 
perturbation mass matrix to a complex symmetric matrix which leads to CP violation. Therefore 
obtain a realistic neutrino mixing matrix with 023 = 45°. We obtain the solar mass splitting, the 
ordering of the neutrino masses is inverted. Using only two sets of the experimental data, we can 
hx all of the parameters of mass matrix and predict the masses of neutrinos and phases. These 
predictions include the following: mi « (4.82 — 4.93)10“^eU, |m 2 | « (4.90 — 5.01)10~^eU, m 3 « 0 
and, (f) ~ (0.687° — 10.31°) as the origin of the Majorana phases. 


1 Introduction 

The results of the neutrino oscillation experiments mm have established that neutrinos have masses. 
The remarkable experimental achievements resulting in experimental data for neutrino give us infor¬ 
mation about neutrino masses and mixing, which can be summarized as follows in Table fTT3ll4]. 


The recent measurement indicate that 613 is non-zero by more than 5cr [ 2 ] and it is small compared to 
the other neutrino mixing angles. 

The lepton mixing matrix in the standard parametrization is given by [S] , 


J 

( C12C13 

S12C13 


/ e^P 

0 


UpMNS = 

— S12C23 — Ci 2 S 23 Sl 3 e®^ 

C12C23 — Si 2 S 23 Sl 3 e®^ 

S23C13 1 

0 

1 

0 

1 

\ S12S23 — Ci 2 C 23 Sl 3 e*^ 

— C12S23 — Sl 2 C 23 Sl 3 e*^ 

C23C13 / 

V 0 

0 
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Parameter 

The experimental data 

The best fit (±lcr) 

Am^JIO-'^eF^] 

7.II-8.18 

7.42 - 7.79 

Amg;^ [10“^eF^] 

2.30-2.65 

2.20-2.54 

2.41 - 2.53 
2.32-2.43 

sin"’ 012 

0.278 - 0.375 

0.307- 0.339 

sin"’ 023 

0.395 - 0.642 

0.405 - 0.639 

0.452- 0.599 
0.535- 0.598 

sin’^ 013 

0.0179-0.0247 

0.0182-0.0250 

0.0201 - 0.0224 

0.0204- 0.0226 

(5 

0 - 27r 

0 

1 

to 


Table 1: The experimental data for the neutrinos’ mixing parameters. When multiple sets of allowed 
ranges are stated, the upper row corresponds to normal hierarchy and the lower row to inverted 
hierarchy. 


where = cosOij and Sij = sinOij (for i,j = (1, 2), (1, 3) and (2, 3)). The phase S is called the Dirac 
phase, analogous to the CKM phase, and the phases p and a are called the Majorana phases and are 
relevant for Majorana neutrinos. 

A Dirac mass term for the neutrinos and charged leptons is written as 

£m = — ^LiM.-D^^VRj + h.C., (2) 

A successful phenomenological neutrino mass model with the flavor symmetry that is suitable for the 
Dirac neutrinos was proposed by Friedberge and Lee (FL) [B]. In this model the mass eigenstates of 
three charged leptons are identified with their flavor eigenstates. Therefore neutrino mixing matrix 
can be simply described by a 3 x 3 unitary matrix U which transforms the neutrino mass eigenstates 
to the flavor eigenstates, {ve, In the pure FL model, one of the neutrino masses is exactly zero, 

partially justifying the smallness of neutrino masses. Morover when p — v symmetry is assumed, the 
matrix U reduces to the experimentally favored Utbm- 

The Dirac neutrino mass operator of the FL model can be written as 

AIfL = a{Vr - {^T -y^i)+b - V^) - Ve) 

+ c{Ve- Vr) {Ve “ Vt) + TUq {VeVe + ■ (3) 

All the parameters in this model {a,b,c and mg) are assume to be real. For toq = 0, this Lagrangian 
has the following symmetry Vg —>■ Vg + z, —>■ Vfg + z, and Vr —>■ Vr + z, where z is an element of 
the Grassman algebra. For constant z, this symmetry is called FL symmetry [ 6 ] in which case the 
kinetic term is also invariant. However the other terms of the electroweak Lagrangian do not have such 
symmetry. The mg term breaks this symmetry explicitly. However we may add that the FL symmetry 
leads to a magic matrix and this property is not spoiled by the mg term. The magic property has 
many manifestations which we shall discus in details. Also it has been reasoned that the FL symmetry 
is the residual symmetry of the neutrino mass matrix after the iS'0(3) x U{ 1 ) flavor symmetry breaking 
[7]. The mass matrix can be displayed by, 

( 6 + C + TOg —b —c \ 

—b a + b + mo —a , (4) 

— c —a a + c + TOg J 

where a oc + Yr^), b oc (Yg^ + Y^g) and c oc {Yrg + Ter) and Yap denote the Yukawa coupling 
constant. The proportionality constant is the expectation value of the Higgs field. It is apparent that 
Mrl possesses exact p — r symmetry only when b = c. Setting h = c and using the hermiticity of Mfl, 
a straightforward diagonalization procedure yields U^^j^^MplUtbm = Diag {mi, m 2 , mg} where, 

mi = 36 + mg m 2 = mo m 3 = 2a+ 6 +mg, (5) 
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and the well-known experimentally favored tribimaximal (TBM) neutrino mixing matrix can be repro¬ 
duced [5]. It has the following form: 


Utbm 



1 

V 3 

1 



V 3 



( 6 ) 


The exact tribimaximal mixing matrix UTBM^ regardless of the model, fixes the element {Ue 3 )TBM = 
0. The values of the mixing angles in the mixing matrix Utbm are consistent with the Table [U except 
013. The role of a non-zero 0i3, or equivalently Ue 3 , is numerus. It is necessary for CP violation in 
neutrino oscillations and may be to explain leptogenesis. For CP violation, of course, both 0i3 and 
the complex phase 6 should be non-zero. Besides, 0i3 ^ 0 will be equivalent to the quark sector that 
mixing between all three generation is a confirmed result, although the mixing angles in the two sectors 
are very different. It is apparent that Eq. m implies the following forms of the neutrino mass matrix 
in the flavor basis, 


TOl 

T 


4 

-2 

-2 \ 

i 

! 1 

1 

1 \ 

1 

( 0 

0 

0 

-2 

I 

1 


1 

1 

1 


0 

I 

-I 

-2 

I 

1 

1 3 \ 

V 1 

1 

1 

1 2 1 

V 0 

-1 

1 


(7) 


where mi (i = 1,2,3) are the neutrino mass eigenvalues. Therefore with respect to the second term 
in Eq. ©, this form of neutrino mass matrix is called democratic matrix m- Democratic matrix is a 
phenomenological Dirac neutrino mass model with S{3)l x S{3)ii flavor symmetry. Democratic basis 
is adopted to produce a flavor-democratic mass matrix in which all matrix elements are equal. In 
this basis, the S 3 operations generate the permutations of three objects, for example, the exchange 
of the first and second indices. The invariance under such transformations require the universal size 
of couplings for three generations if they belong to a three-dimensional representation of S 3 and the 
Higgs held is in the singlet. 

The smallness of 0 i 3 compered to the other mixing angles persuade us to modify the neutrino 
mixing matrix by a small perturbation on the basic tribimaximal structure and could lead to a realistic 
neutrino mixing matrix. We focus our interest on the neutrino mass model was proposed by Friedberg 
and Lee (FL) in a very special case which leads us to obtain the Democratic mass matrix in the havor 
basis. Democratic matrix is suitable to generate 0 i 3 7 ^ 0 by starting from an initial tribimaximal 
form. However at the hrst sight the resulting of its mass spectrum with the Utbm is experimentally 
unfavored. Therefore the purpose of this paper is to hnd the experimentally favored mass spectrum of 
the Democratic matrix along with a realistic neutrino mixing matrix by a small perturbation on the 
basic tribimaximal structure. There are lots of neutrino mass models from which one can obtain the 
tribimaximal form of the neutrino mixing matrix . As well many theoretical and phenomenological 
works have discussed massive neutrino models to generate 0 i 3 0 in different ways starting from an 

initial tribimaximal form m- 

In order to have CP-violation in the standard parametrization given in Eq. o, the necessary 
condition is (5 7 ^ 0 and 0 i 3 ^ 0. There are four independent CP-even quadratic invariants, which can 
conveniently be chosen as Du, 1 /^ 31 / 13 , C/^i 1/21 and U 23 U 23 and three independent CP-odd quartic 
invariants m, 


J = Im{UnU* 2 U*^U 22 ) 

h = Tm[{Ul^Ui2f] 

h = Xm[{Ut^Ui3?]- ( 8 ) 

The Jarlskog rephasing invariant parameter J m, is relevant for CP violation in lepton number 
conserving processes like neutrino oscillations. Ii and I 2 are relevant for CP violation in lepton number 
violating processes like neutrinoless double beta decay. Oscillation experiments cannot distinguish 
between the Dirac and Majorana neutrinos. The detection of neutrinoless double beta decay would 
provide direct evidence of lepton number non-conservation and the Majorana nature of neutrinos. 
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Many theoretical and phenomenological works have discussed massive neutrino models that break 
/r — T symmetry as a prelude to CP violation[TB]. 

In this paper we generalize the FL model by introducing complex parameters which can ultimately 
be linked to complex Yukawa coupling constants. We concentrate on the massive FL Dirac model, 
imposing the obvious constraint that mass eigenvalues be real. In our model there is a restricted 
region of the parameter space where CP violation occurs. We focus on a special point on the border of 
this region which create democratic neutrino mass matrix in the parameter space. Democratic neutrino 
mass matrix is experimentally unfavored, due to mi = m 3 = 0. Therefore, we consider a case with 
mi 0 by adding the symmetry breaking term based on the FL symmetry. Afterward on the basic 
TBM structure, the universality of the experimental data insert some constraints on the elements of 
the mass matrix. Therefore we obtain our unperturbed neutrino mass matrix with /r — r symmetry and 
magic property in the flavor basis. We find that the solar mass splitting is absent, hence we generate 
this splitting via a small perturbation which is responsible for Ue 3 ^ 0. This work can be broken down 
to two clearly distinguishable stages. In the first stage, we obtain the perturbation mass matrix by 
the help of the perturbation theory in the mass basis with real parameters which breaks mildly the 
/i — T symmetry and this may generate ^13 ^ 0 but will have no CP violation. In the second stage, 
we extend our work to the case of CP violation, therefore complex perturbed mass matrix have been 
generated. So we obtain nonzero values for both ^13 and S, and this leads to CP violation along with 
a realistic neutrino mixing matrix. 

The main outline of the paper is as follows. In section 2, we introduce our model which obtain 
from a special point in the generalized FL model and show the results of the aforementioned stages, 
separately. We find our symmetry breaking and perturbation mass matrix in the flavor basis which 
causes CP violation and solar mass splitting, moreover we get realistic neutrino mixing matrix. In 
section 3, we compere two sets of the experimental data onto the allowed region of our parameter 
space. We then find the allowed region for the perturbation parameters. Also we predict the masses 
of the neutrinos and phases. We then not only check the consistency of all of our results with the 
available experimental data, but also present our predictions for the actual masses. In section 4, we 
state our conclusions. 


2 The Model 


In this section, we generalize the FL model by adding complex Yukawa coupling constants, in order 
to obtain CP violation. This is accomplished by obtaining non-zero values for sindi 3 and S. However, 
we demand the eigenvalues of the mass matrix to be real. We find that only one particular choice 
allows for minimal breaking oi fi — t symmetry, i.e. (a G 5ft; 6 , c £ C and b = c*). This requirement 
leads to a non-hermitian mass matrix. For simplicity of notation we define the parameters as follows: 
5ft ( 6 ) = 5ft (c) = br and S (&) = -S (c) = 43- 
The neutrino mass matrix M' is given by 


M' = 


2br -\- Too 

— br 

— br 

—br a 

+ br + mo 

—a 

— br 

—a 

a + br + mo 

( 0 -1 

1 \ 


-1 1 

0 . 


1 1 0 

-1 / 



(9) 


Notice that M[, and are both magic and symmetric matrices since one of the eigenstates is 

and we choose it to be \i^ 2 ) in order to be consistent with Eq. ([5]). 

^ The parameters indicating the measure of CP violation and fi — t symmetry breaking turn out to be proportional 
to B and therefore we expect it to be small. 
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A naive diagonalization of M'^ yields, 


rhi 

= (a -I- 2br + mg) + ^/(a — br)^ — 3B^ 


rfi2 

= mg. 


m 3 

= (a -I- 2br + mg) — \/(a — 5^)^ — 35^. 

( 10 ) 


However the usual diagonalization procedure is correct only for hermitian matrices, where a simi¬ 
larity transformation by a unitary operator, i.e. diagonalizes the matrix. Therefore, 

the results indicated in Eq. (uni) are correct only in the limit H —>■ 0 , and the results in this limit suffice 
for our analysis to follow. Comparing our results in this limit with those shown in Eq. ([5]) , we conclude 
a <b. 

Since M' is a non-hermitian matrix, we need two distinct unitary matrices U and V to diagonalize 
itl we do not display the explicit form of U and V and only mention that V = U*. The resulting 
correct diagonal matrix is obtained by = WM'^V and its elements are as follows, 


m'l 

m'2 

m 3 


iB{a 


mo, 

iB{a 


6 r) “t“ -t- (n -t- 267 . -f mo)^ — (a — bj- -t- (a -t- 2 .br -I- mg)^ 

a + 2{br + iB) -\- mo 


&r) “f 3B^ (tt -t- 2bj- mg)^ “t“ (a — br iB') ^3B‘^ -t- (n -t- 2br mg)^ 
a + 2{br + iB) + mg 


( 11 ) 


After the diagonalization we find that only m'l and m 3 are complex. We can extract the phases 
and transfer them to the mass eigenstates in the Dirac case [14]. The most general form the diagonal 
mass matrix can be written as, 




( 12 ) 


In our model a automatically turns out to be zero. We would dispense with the overall phase even 
if it was not zero. Using the fact that m^ is real, we obtain /3 = 7 . This implies that the aig{m[) = 
— arg(mg) = 2)^ Using these conditions in Eq.(ll) we obtain 


B = ± 


— (2a + br + mg) (36^ + mo) 


(13) 


Erom the requirement of the reality of B we obtain —^<br< —(2a + mg). Notice that the lower 
bound of br is simply a check on the condition mi > 0. The requirement that in the limit U —>■ 0, U 
and V should both approach Utbm given in Eq. ([5]), yields 2br -I- a-I- mg > 0. Combining the condition 
for reality of B with a < b, we obtain 3a -I- mg < 0. From this and the overall symmetry of the F.L 
model we conclude that br < ( 0 . Notice that the occurrence of CP violation is possible only in a 
restricted region in the a-br plane where B ^ Q. Figure o illustrates the region of the parameter 
space where CP violation occurs. 


We have examined all details about CP violation region, obtainable consequences of this region and 
consistency with the experimental data in Ref |15j . 

In Figure (|T|) the right angle of the dark triangle is a special poin that display a = br = —^ at 
which B = 0. Therefore in our model we cannot have CP violation at this point. In this paper, we 

^ These matrices can be easily obtained by diagonalizing and Ml, separately. U and V are the conven¬ 

tional transformation matrices for the left-handed and right-handed neutrinos, respectively. 

®Note that disregarding the overall phase amounts to the following: is real and Det{U) = iQil. 

This conclusion is consistent with the results of experiments on solar neutrino oscillation which indicate that 
m2 > mi. 
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Figure 1: CP violation is possible only in the right-angled triangle in our parameter space. The axis are 
defined by q = and The dark triangle displays the allowed region within our model. (The line 

above the base of the triangle is given by 2br -I- a -b mo = 0) 

focus our interest only on this point, where the neutrino mass matrix reduces to the following special 
form: 


M' 


/ 1 0 0 

(36r + mo) I 0 1 0 

\ 0 0 1 

(3br + mo)l — brD. 


— br 


1 1 1 
1 1 1 
1 1 1 


(14) 


the first term in Eq. (IT41) is identically zero and reduces to —brT>. Therefore M^, gets the Democratic 
form that can be realized by S 3 permutation family symmetry [TU]. Here we separately impose the S 3 
permutation symmetry for the left- and right- handed neutrinos as a flavor symmetry. They are shown 
by S 3 L and Son, respectively. Therefore leads to invariant M' under S 3 L x S 3 R symmetry. 

M[, can be diagonalized by the Utbm matrix which yields m 2 = mo and mi = m 3 = 0. Therefore 
the model is experimentally unfavored. The possible way to rescue the model is to obtain the non¬ 
vanishing mi or m 3 . We consider a case to have non-zero mi by adding the symmetry breaking term 
based on the FL symmetry in what follows. 

The combination of the FL symmetry with the fJ- — t symmetry is a kind of translational family 
symmetry that is called twisted FL symmetry [IT]- The twisted FL symmetry for the Dirac neutrinos 
impose on the left-and right-handed neutrinos separately as follows, 

l^Li VLi ~ ^ij^Lj + ^LjZ 

^Ri = + ^Ro^ > ( 15 ) 

where z is a space-time independent Grassmann parameter, = 0, A = (Ai, A 2 , A 3 )^ are c-numbers, 
and S is the permutation matrix between the second and third families: 

1 0 0 \ 

0 0 1 . ( 16 ) 

0 10 / 

Two independent /i — t permutation symmetries make the Dirac mass matrix as 

/ D -2C -2C \ 

= -2B -A -A \ ■ 

\ -2B -A -A ) 

while the translational symmetries in Eq. (Ha lead to the conditions Xrj = 0 
The resulting form of the mass matrix depends on the correlations among Al^ and 
we have assumed the uniform translation, that is r]Li oc (1, 1 ,1) and rijn oc ( 1 ,1, 1 ). 


(17) 

and Aii M/’ = 0 . 
Xm- In this work, 
Then, under these 
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transformations, the breaking mass terms for the S^l x flavor symmetry of the Dirac neutrino 
takes the form, 


(18) 


We add the mass matrix, M^, in Eq. (IT^ to Eq. (fHl) as the breaking term for the S 3 L x S 3 R flavor 
symmetry. Consequently, the total Dirac neutrino mass matrix is written as 



mbd ^ 


/ 1 

1 

^ \ / 

' 4 

-2 

-2 

1 

1 

1+5 

-2 

1 

1 

\ 1 

1 

1 / \ 

-2 

1 

1 


(19) 


We call It Broken Democtratic (BD) neutrino mass matrix and Its admissible neutrino mass spectrum 
is given by 


fhi = 6 g, m 2 = mo, m 3 = 0 , 


( 20 ) 


which suggests the inverted neutrino mass hierarchy and, the TBM mixing matrix can be also realized. 
This is a minimal scheme: the twisted EL symmetry for the Dirac neutrinos is introduced as the 
symmetry of the breaking term. We note that the whole mass matrix in Eq. (HI violates the S^l x S^n 
symmetry but still preserves the fJ. — t symmetry. 

Now on, we want to find a small perturbation mass matrix along with a realistic neutrino mixing 
matrix from BD mass matrix on the basic TBM mixing matrix. Therefore from Eq. ([7]) a general mass 
matrix, satisfying Utbm when represented In the flavor basis and the most general form of it is. 


M.V = Utbm 


A 



B A+C B-C 
B B-C A+C 


TBM 


( 21 ) 


and we have set 


m = 


T,rni 


A; 


3j 


= (m3 - ruj), for j = 1 , 2 . 


Therefore by Eq. (I7T|) and Eq. (1771) , we obtain A = m — , B — 


A31 —As 


and C = 


( 22 ) 

Substituting 


the m^s in Eq. (ITHl) into the expressions in Eq. (1771) . we obtain m = ^+2g, A 32 = —mo and A 31 = — 6 g. 
We work in a flavor basis in which mixing in the lepton sector is determined entirely by the neutrino 
mass matrix. 

The experimental data have now dehnltely confirmed that Am^i = rn\ — m\ > 0 and too tiny. By 
substituting m^s in Eq. (I7D1) into Arri^i, we obtain that Then we approximate A 32 — A 31 = A 

thus A ~ —mo ~ — 6 ( 7 , which indicate that 5 ~ Since A is negativ^l. 

BD neutrino mass matrix with these approximations in the flavor basis, the unperturbed neutrino 
mass matrix, is given by 

m — y 0 0 

M,? ~ I 0 
0 


m + ■= —■= 


By substituting the expression g ~ 
(0) (0) 


mn 

6 


m + 

f(0) 


A 


(23) 


ml ~ = m —— = mo 


the eigenvalues of Mj) ' are, 

A . (o)_ 2A 


and mq = m -\ -~ 0 . 

3 


( 24 ) 


^ A is important quantity because the scale for atmospheric neutrino oscillations is set by A. 
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It is apparent that in Eq. (1^ . as before, possesses exact ^ — t symmetry and the magic property. 
Therefore, the mixing matrix is still Utbm- ^i, and m 2 are real and positive, although at this level, 
the solar mass splitting is absent. As before m 3 = 0, therefore the ordering of neutrino masses is still 
inverted. Our aim is to obtain the solar mass splitting by the same mass perturbation that is cause 
of 013 7 ^ 0 and CP violation. Finally, in this model CP violation conditions necessarily mandate that 
fj. — T symmetry should be broken. An interesting question is whether 023 = 45° holds after the fJ- — t 
symmetry breaking. 

Our method for finding the structure of the neutrino mass matrix, is to use perturbation theory. 
So that Ml, = + Mjf , where << M°. Both and will be symmetric and could be, 

in general, be complex. However, as obtained in Eq. (1231) from BD neutrino mass matrix and the 
tribimaximal mixing feature is real and symmetric, therefore it is Hermitian. We will consider the 
cases of real and complex M,f in different stages. Therefore the result of our work, the solar mass 
splitting, 013 7 ^ 0 and, CP violation will be obtained. 

The eigenstates of in the mass basis, the unperturbed mass eigenstates, are found useful. These 
are as follows: 


= ( [J j ’ " ( 0 ) ’ " ( 1 ) ’ 

which the first two mass eigenstates are degenerate. We investigate such that and 1 ^ 2 °^ 
its nondegenerate eigenstates: |M,f = m[^'^Sij where {i,j = 1,2), with 7 ^ Then 

in this mass basis, we take {M ^)33 = 0 and we only consider {M^)i 3 and {M^) 23 . Therefore the 
basis vectors Vi and 1^2 are chosen with the aim that they reproduce the correct solar mixing and the 
physical basis is fixed by the perturbation. Needless to say, that when eigenstates in Eq. (1251) expressed 
in the flavor basis are the columns of Utbm in Eq. ([HI). Therefore in the flavor basis eigenstates are as 
follows: 


'1">) = 


2 
3 

z 

V Ve 


Z) = 


v4 


,( 0 ) 


) = 


0 

1 

C2 


(26) 


At the first stage our aim is to obtain the perturbed mass eigenstates that when expressed in 
the flavor basis are the columns of the matrix in Eq. ©• For this, at first, we choose which is 
symmetric, to be real and therefore Hermitian. Hence, this may generate a nonzero 0 i 3 but will not 
lead to CP violation since this computation necessarily yields (5 = 0. For the perturbation expansion 
we keep terms up to linear in S 13 . To first order we have 


3^3 


(27) 


Where, 

C37 = -Cj3 = (to3°^ - mf'')-^< >, {j ^ 3 ). ( 28 ) 

In this case, the coefficients Csj are real and in the mass basis is proportional to p7] . 

\v 3 ) given in Eq. (I?7ll should correspond to the third column of the mixing matrix Upmns in Eq. ([T]) 
with (5 = 0 . By using Eq. (Eli) in the flavor basis C31 and C32 are easily determined. Therefore explicitly 
we get the matrix equation, as follows: 






Readily one obtains, to order linear in S 13 , C 31 = — \/|si 3 and C 32 = \ 3 S 13 , where maximality of the 


2-3 mixing angle has been used (023 = 45°). Therefore in the mass basis by using Eq. (IMl) and Eq. ( 
could write, )^3 = and = —moy^sis. 

We continue this discussion in the second stage as follows in order to obtain CP violation. We let 
be complex symmetric matrix but not Hermitian and it is the same for the total mass matrix 
Ml , = -I- . Therefore, this is accomplished by obtaining nonzero values for sin 0 i 3 and (5^0. 

The columns of the mixing matrix U in the Eq. (U]) are eigenvectors of + 

where we have dropped a term, which is To proceed, we recall that unperturbed 

is Hermitian and its eigenstates are the same as the columns that produce Utbm, in the 
Eq. dni), and its eigenvalues are |m 2 °^p, and In place of Eq. (l2^ we have 


C3, = -C*. = 


J3 




) Mj3, (j ^ 3) 


(30) 


where Mjs = < > and \v 3 ) be reproduced to first order by substituting 

the expression for Cgj given by Eq. (1501) into the Eq. (|07)) . Therefore by using the appropriate variant 

of Eq. (l29)) . we get in this case C31 = — and C32 = Notice that by using the 

symmetric nature of in the mass basis simply could relate the C 31 and C 32 to the elements of 
perturbation as follows: 






(0) 


> 


= mf\M^)j 3 + {j ^ 3) 


( 0 )/ 


(31) 


By using Eq. (|5T|) could write (M^f )i 3 = toq 
the expression for m‘f\ given by Eq. (IMl) . can find 


|si 3 e*'' and (M')23 = -mo 


3 S 13 e*’* that by using 


77 =-tan WtanJ), and S13 =/(??)si3, 


(32) 


where 


fijl) = 


(to^°^)^ -|- Y + 2to^°^TO3°^ cos2?7 


1/2 


,(o) 


.( 0 ) 


(33) 


The formulae in Eq. (1321) display that S 13 < S 13 , in our model = 0, therefore the equality 
holding. The range of ry and S is the same as each other which is {0, 27r}. 

Therefore we have focused on finding 0 i 3 ^ 0 through a perturbation starting from the tribimaximal 
neutrino mixing matrix. Now we investigate the solar mass splitting. A straightforward calculation in 
our model yields (M^)i 2 = {M^) 2 i = 0. The first order corrections to the neutrino mass are obtained 
from 5ij = < 

this order: 

= 0 and m^^^ ^ 0. (34) 

Therefore in the mass basis this implies that {M ^)22 7 ^ 0 and other diagonal elements of the 
perturbation matrix are zero. Such a correction display a nonzero solar mass splitting, in which 


vy >. We seek that the first-order of neutrino mass corrections arise at 




' = 1712 — mi, and Am^i = (m 2 )^ — (mi)^ is positive. 


Therefore we write the hnal perturbation matrix in the mass basis, as follows: 


/ 


MT = - 


mo Si3 


0 
0 




e 


Ipir] 

3^ 


(35) 
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and 



Too si3 ' 


( 36 ) 


The solar mass splitting fixes the dimensionless parameter e. In general the solar mass splitting can be 
complex. It mentioning that the Majorana mass 7712^^ = exp(i(p). If we write TO 2 = +m 2 ^^ = 

|to 2 | exp(i<^) then one has, 


|mi| 

|m2| 


= Too, Imal = = 0, 

+ (|to 2^^|)^ + 2777-2°^|rn2^^I cos(^ 


= tan 


-1 


m 


w. 


sm(/? 


.(0) 


.(i)i 


COS(f 


■ (37) 


For the Dirac neutrinos, these phase (j) can be removed, and for the Majorana neutrinos, these 
phase remain as Majorana phase and contribute to CP violation. 

Using degenerate perturbation theory [28], to order linear in S 13 , and in Eq. (I35|) we obtain 
the neutrino mixing matrix with <5 ^ 0 , as follows: 



The nonzero 5 shows CP violation in the lepton sector. The same form of U in Eq. (1551) has been 
discussed in |29j from a different motivation and its consistency with the observed mixing angles noted. 

To relate perturbation mass matrix to the perturbed mass matrix, BD, we rewrite in the flavor 
basis. Therefore, by using Eq. we obtain perturbation mass matrix in the flavor basis as follows: 


Muf = -iTT'O Sl3 



0 

-1 

1 


-1 1 

0 0 

0 0 


£ 



1 1 
1 1 
1 1 



(39) 


The first term on the right-hand side is responsible for ^13 and the second for Arri^i- 

A rephasing-invariant measure of CP violation in neutrino oscillation is the universal parameter J 
m given in Eq. (|5|), and it has a form which is independent of the choice of the Dirac or Majorana 
neutrinos. Using Eq. (1551) and Eq. (1551) the expression for J simplifies to. 


J = = 

Notice that both S 13 and ry have to be non-vanishing in order for CP violation to be present in the 
lepton sector. 


3 V 2 


0)13 '/• 




3 Comparison with experimental data 

In this section we compare the experimental data with the results obtained from our model. Therefore 
we obtain the allowed ranges of the neutrino masses along with the perturbation term. We do this in 
two stages. In the first stage, we compere the results of neutrino masses in our model in Eq. (1371) with 
the mass constraints of which obtained from the experimental data, in Table [TJ Therefore the 

allowed ranges for mi = toq is obtained as follow, mi = toq ~ (4.82 — 4.93)10“^eU, which agree well 
with the allowed ranges for mi and mg that obtained in Ref [la¬ 
in the second stage, we expect that the different nonzero terms of the perturbation matrix Eq. (1351) 
are roughly of similar order. We may then expect £ ^ C>(1), and could predict the order of m^'^ 
in Eq. (1551) . Therefore by using mg 0(10“^) from previous stage and sin 6(13 ^ 0(10“^) from the 
experimental data in Tabled] we predict m^'' ^ 0(10“^). As a result of this prediction we obtain 
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the allowed ranges of the neutrino perturbation parameters, and (p. Therefore the allowed range 
of m 2 will be obtained, too. We do this by mapping neutrino mass constraints obtained from the 
experimental data, Am^i, onto our parameter space. In Figure (I2|) we have plotted the overlap regions 
of and our model along with the allowed ranges of mo onto perturbation parameter space, \m^'^ \, 

and (/?. 



Figure 2: In this figure the whole region of the plane which is allowed by our model along with the 

allowed ranges for mo is shown. The overlap region of the experimental values for Am^i. In the zoomed box 
we have magnified the overlap regions. 


The results of the perturbation parameters are as follows. 


di)| 




(0.001 - 0.009)ey, 
±(36° -90°). 


(41) 


In this work, we predict the neutrino masses and the phases. The results are as follows, which all 
of them agree well with the allowed ranges of the neutrino masses that obtained in Ref [15] . 
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( 42 ) 


TOi « (4.82-4.93)10"^ey, 

Imzl ~ (4.90-5.01)10-^ey and (/> Ri (0.687° - 10.31°), 

TO3 Ri 0. 

As mention before (j) is the origin of the Majorana phases which arise from the perturbation. We 
would dispense with the overall phase, exjp^icj)), therefore we have two phases that appear in the mass 
eigenvalues shown in Eq. (H^ . i.e., p = (z^) and a = exp[—z(^ — For the Dirac neutrinos, these 
phases can be removed, and for the Majorana neutrinos, these phases remain as Majorana phases and 
contribute to CP violation. 

An important experimental result for the sum of the three light neutrino masses has just been 
reported by the Planck measurements of the cosmic microwave background (CMB) at 95% CL |34) . 
which is 


y^m^< 


0.230 eV (PlancA:+WP+highL+BAO) 
0.247 eV (F/ancA:+WP+BAO) 

0.663 eV (FZancA:+WP+highL) 

0.933 eV (Planck+WP) 


(43) 


This sum in our model is k, (0.097 — 0.099)eE, which is exactly consistent with the constraint 

of the above. 


4 Conclusions 

In this paper we have proposed a generalization of Friedberg-Lee neutrino mass model, in which CP 
violation is possible. In our model the coefficients are allowed to be complex, with the constraint that 
the mass eigenvalues be real. We find and display the region in our parameter space where CP violation 
is possible. However, we have focused only on a special point on the border of this region which create 
a neutrino mass matrix with permutation family symmetry. Democratic matrix, in the parameter 
space. On the basic Utbm, Democratic matrix is an experimentally unfavored due to its degenerate 
mass spectrum. Therefore, we improve our model to obtain the experimentally favored neutrino mass 
spectrum by the breaking mass term which preserve the twisted FL symmetry. Therefore, we obtain 
a mass spectrum without degeneracy, the model predicts the massless third generation neutrino and 
inverted mass hierarchy. The tribimaximal mixing matrix can be also realized. We obtain a realistic 
neutrino mixing matrix along with a small perturbation mass matrix from the BD mass matrix on 
the basic Utbm- We start with the BD mass matrix that the experimental data are caused some 
constraints on it. Therefore the resulting mass matrix is symmetric with p — t symmetry, and has 
magic feature, however it loses the solar neutrino mass splitting. At this level, our study consists of 
two stages; in the first stage, we obtain the elements of the perturbed mass matrix in no CP violation 
case by the help of perturbation theory. In the second stage, we extend our study to the case of CP 
violation, therefore obtain the complex elements of the perturbed mass matrix in the flavor and mass 
basis which break softly the p — t symmetry and magic feature. We get a realistic mixing matrix with 

5 0 and 023 = 45°. Also obtain solar neutrino mass splitting. Comparing the results of our model 
with Am'^i from the experimental data, we obtain the allowed range of mi « (4.82 — 4.93)10“^eF. 
Also we predict the order of solar mass splitting is ^ 10“^. Mapping the allowed ranges of Am^i from 
the experimental data, along with the allowed range of mo onto perturbation parameter space almost 
pinpoints the values of our parameters. Then we predict |m 2 | ~ (4.90 — 5.01)10“^eF, as well inverted 
hierarchy is preserved with m 3 = 0. Our predictions for the neutrino masses are to tested by future 
experiments. 

®Note that disregarding the overall phase amounts to the following: Det{U) = 1 I3(JI . li (f) = ^. 
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